calculate_hum_stats <- function(row, temp_data) {
  # Adjust temp_data with data.table() format
  setDT(temp_data)
  
  # Dates are dates 
  if (!inherits(row$date_start_week, "Date")) {
    row$date_start_week <- as.Date(row$date_start_week)
  }
  if (!inherits(row$date_end_week, "Date")) {
    row$date_end_week <- as.Date(row$date_end_week)
  }
  
  # Extract values for each row
  start_date <- row$date_start_week
  end_date <- row$date_end_week
  comuna_id <- row$comuna
  
  # New temporal data with filter by time and district
  week_temperatures <- temp_data[fecha >= start_date & fecha < end_date & comuna == comuna_id, 
                                 .(h_min, h_max, xh_diario)]
  
  # Blank Results 
  result <- as.list(row)
  
  # NA with missing values when haven't information in temp data 
  if (nrow(week_temperatures) == 0) {
    result <- c(result, list(hum_mean=NA, hum_min=NA, hum_max=NA,
                             hum_mean_min=NA, hum_mean_max=NA))
    return(as.data.table(result))
  }
  
  # Calculate stats and return results as data.frame()
  result <- c(result, list(
    hum_mean = mean(week_temperatures$xh_diario, na.rm = TRUE),
    hum_min = min(week_temperatures$h_min, na.rm = TRUE),
    hum_max = max(week_temperatures$h_max, na.rm = TRUE),
    hum_mean_min = mean(week_temperatures$h_min, na.rm = TRUE),
    hum_mean_max = mean(week_temperatures$h_max, na.rm = TRUE)
  ))
  
  return(as.data.table(result))
}

calculate_temperature_stats <- function(row, temp_data) {
  # Asegurarse que tanto row como temp_data son data.table
  setDT(row)
  setDT(temp_data)
  
  # Convertir fechas solo si es necesario
  row[, date_start_week := as.Date(date_start_week)]
  row[, date_end_week := as.Date(date_end_week)]
  
  # Filtrar los datos de temperatura
  week_temperatures <- temp_data[fecha >= row$date_start_week[1] & fecha < row$date_end_week[1] & comuna == row$comuna[1]]
  
  # Resultados en blanco o con NA si no hay datos
  if (nrow(week_temperatures) == 0) {
    return(data.table(row, temp_mean=NA, temp_min=NA, temp_max=NA, temp_mean_min=NA, temp_mean_max=NA))
  } else {
    return(data.table(row, 
                      temp_mean = mean(week_temperatures$xt_diario, na.rm = TRUE),
                      temp_min = min(week_temperatures$t_min, na.rm = TRUE),
                      temp_max = max(week_temperatures$t_max, na.rm = TRUE),
                      temp_mean_min = mean(week_temperatures$t_min, na.rm = TRUE),
                      temp_mean_max = mean(week_temperatures$t_max, na.rm = TRUE)))
  }
}

glimpse(births)

length(unique(births$id))

cat("% de Avance Muestra:", round((length(unique(births$id))/361922)*100,0), "%")

births <- births %>% 
  mutate(temp_mean=if_else(is.infinite(temp_mean), NA, temp_mean),
         temp_min=if_else(is.infinite(temp_min), NA, temp_min),
         temp_max=if_else(is.infinite(temp_max), NA, temp_max)
         )

temp <- births %>% 
  group_by(zona, ano_nac, mes_nac, dia_nac) %>% 
  summarise(temp_media = mean(temp_mean, na.rm = TRUE), 
            temp_min = mean(temp_min, na.rm = TRUE), 
            temp_max = mean(temp_max, na.rm = TRUE), 
            peso = mean(peso, na.rm = TRUE)) %>% 
  mutate(date_nac=make_date(year=ano_nac, month=mes_nac, day=dia_nac)) %>% 
  pivot_longer(!c(zona, ano_nac, mes_nac, dia_nac, date_nac, peso), 
               names_to = "temperatures", 
               values_to = "values")

peso <- births %>% 
  group_by(zona, ano_nac, mes_nac, dia_nac) %>% 
  summarise(peso = mean(peso, na.rm = TRUE)) %>% 
  mutate(date_nac=make_date(year=ano_nac, month=mes_nac, day=dia_nac))
  
ggplot(temp, aes(x = date_nac)) +
  geom_line(aes(y = values, fill=factor(temperatures), color=factor(temperatures))) + 
  scale_y_continuous(breaks=seq(-10, 40, by=10)) +
  # geom_line(data=peso, aes(y = peso/100)) + 
  # scale_y_continuous(name = "Temperatura", limits = c(-20 , 60),
  #                    sec.axis = sec_axis(~.*100, name = "Peso del nacimiento (kg)")) +
  scale_x_date(date_breaks = "1 year", date_labels = "%Y") +
  facet_wrap(~ zona, scales = "free_x", ncol=1)  +
  labs(x = "", y = "Temperatura") +
  theme_light() +
  theme(axis.text.x = element_text(angle = 0, hjust = 1), 
        legend.title = element_blank(),
        legend.position = "top", 
        strip.background = element_blank(),
        strip.text = element_text(color="black", face = "bold"))